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ABSTRACT 

Context. Asteroseismology has great potential for the study of metal-poor stars due to its sensitivity to determine stellar ages. Solid 

detections of oscillation frequencies in stars with well constrained fundamental parameters, combined with a known rotation period, 

should significantly advance our understanding of stellar structure and evolution in context with metallicity effects. 

Aims. Our goal was to detect p-mode oscillations in the metal-poor sub-dwarf 85 Peg A and to search for other variability on longer 

timescales. 

Methods. We have obtained continuous high-precision optical photometry of the binary system 85 Pegasi with the MOST 
(Microvariability & Oscillations of STars) space telescope in two seasons (2005 & 2007). The light curves were analyzed using 
traditional Fourier techniques. Furthermore, we redetermined v sin i for 85 Peg A using high resolution spectra obtained through the 
ESO archive, and used photometric spot modeling to interpret long periodic variations. 

Results. Our frequency analysis yields no convincing evidence for p-modes significantly above a noise level of 4 ppm. Using sim- 
ulated p-mode patterns we provide upper RMS amplitude limits for 85 Peg A. After removal of instrumental trends the light curve 
shows evidence for variability with a period of about lid and this periodicity is also seen in the follow up run in 2007; however, 
as different methods to remove instrumental trends in the 2005 run yield vastly different results, the exact shape and periodicity of 
the 2005 variability remain uncertain. Our re-determined v sin i value for 85 Peg A is comparable to previous studies and we provide 
realistic uncertainties for this parameter. Using these values in combination with simple photometric spot models we are able to re- 
construct the observed variations. 

Conclusions. The null-detection of p-modes in 85 Peg A is consistent with theoretical values for pulsation amplitudes in this star. The 
detected long-periodic variation in the 85 Peg system must await confirmation by further observations with similar or better precision 
and long-term stability. If the lid periodicity is real, rotational modulation of surface features on one of the components is the most 
likely explanation. 

Key words, techniques: photometric - stars: individual: HD 224930 - stars: individual: 85 Peg - stars: oscillations - stars: rotation - 
starspots - binaries:general 



1. Introduction 



With observations dating back to the mid-19 th century by Baily 
(1 18451) . 85 Pegasi (HD 224930) is a very well studied visual bi- 
nary system with an angular separation of 0.83 arcseconds. Due 
to its proximity (d = 12 pc) and brightness (V = 5.75), the or- 
bital elements have been det ermined to hig h accuracy through 
photometry and spectroscopy. iGriffinl (12004) provides a compre- 
hensive review of the wealth of data obtained for the 85 Peg sys- 



Send offprint requests to: Daniel Huber, huber@astro.univie.ac.at 

* Based on data from the MOST satellite, a Canadian Space Agency 
mission, jointly operated by Dynacon Inc., the University of Toronto 
Institute for Aerospace Studies and the University of British Columbia, 
with the assistance of the University of Vienna. 



tem, including a spectroscopic orbital solution with an orbital 
period of 26.31 years and a total mass of 1.49 Mq. The dynam- 
ically measured relative masses of the binary system, however, 
indicating a secondary equally or only slightly less massive than 
85 Peg A, disagree with the observed magnitude difference of 
Amy 3 between the components. It was suggested as early 
as 1948 bv lHalll d 19481) that the secondary star i n 85 Peg may be 
itself a binary system. Evolut iona ry models by iFernandes et ail 
d2002l) . lD , Antona et all d2005l) and lBach et ail d2007h all support 
the binarity of 85 Peg B. 

The primary 85 Peg A is metal-poor, with mea surements 
of fFe/HI = -0 88 dAllende Prieto et alj 120041) . -0.78 
(iHolmberg et alj2007l) and [M/H] = -0.69 (van't Veer 2000; cf. 
IFernandes et alj2002l) . The space velocity of the star, along with 
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Table 1. Main parameters of 85 Peg A and B compiled 
from the literature. Val ues in brackets for c omponent B 
are suggested values by ID' Antona et al.l (120051) for the sus- 
pecte d binary system 85 P egBa and Bb. Reference co ding: 
m Iten Brummelaar et al. (2000) (2) [Fulbrightl J2000h. (3 ) 
Gri ffinl (12004. (4)lD' Antona et al](l2005l) . (5) lBach et all (120071) . 
(6) iFernandes et al l (l2002h. (7) van't Vee r et a l. (2005); cf. 
ID' Antona et all (120051) . (8) lHolmberg etal] (120071) 



parameter 


85 Peg A 


85 PegB [a,b] 


Ref 


V [mag] 


5.81+0.03 


8.89 + 0.29 


1 


Sp.Type 


G5 IV 


K6-8V 


2 


M/M Q 


0.77 ±0.05 


0.72 + 0.05 [0.5,0.2] 


3,4,5 


L/L e 


0.617 + 0.02 


0.072 + 0.03 [0.064,0.0085] 


4,6 


7W [K] 


5600 ± 50 


4200 ± 200 


6 




4.6 + 0.1 


4.8 + 0.2 


4,7 


[Fe/H] 


-0.78 




8 


age [Gyr] 


9.3 + 0.5 
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its low metallicity, qualify i t as a sub-dwarf. 85 Peg A was clas- 
sified as a G5 sub-dwarf bv lFuihrig ht (2000). Based on only the 
visual magnitude difference (neglecting the contradicting mass 
ratio from the orbital solution), 85 Peg B is usually assumed to 
be a K6-8 dwarf. 

IFernandes et al.1 (|2002|) modeled the initial abundances, age 
and mixing-length parameters of 85 Peg A and estimated an age 
of 9.3 +0.5 Gyr. Using effective temperatures der ived from spec- 
trosco py and the measured Hipparcos parallax, ID' Antona et all 
(2005) computed evolutionary tracks for a set of possible metal- 
licities and ages of the system. Their evolutionary models al- 
low ages between about 8 and 14 Gyr. They generated pulsa- 
tional eigenspectra from their models, predicting that solar-like 
p-modes of low degree {{ = 0-3) in 85 Peg A could occur across 
a frequency range of 1-5 mHz (^ 80^1-OOd 1 ). Table [TJ lists 
all important parameters concerning the binary system 85 Peg 
which are relevant for this paper. 

The detection of oscillations in stars like 85 Peg A is par- 
ticularly interesting due to the sensitivity of asteroseismology 
to determine stellar ages. While the small frequency separa- 
tion is the more sensitive observable for this purpose, even a 
measurement of the large separation considerably constrains 
the parameter space which often is more uncertain for metal- 
poor objects than for other stars. The first successful applica- 
tion of asteroseismology t o constrain the age o f a metal-poor 
star has been performed by Beddi ng" et al.1 (|2006) using Doppler 
velocity observatio ns of the subgiant vlndi ([Fe/H] = -1.4). 
Even before that, Mos ser et al.1 Q2005) detected oscillations in 
velocity measurements of HD 49933 ([Fe/H] = -0.4) which 
were confirmed photometrically by the CoR oT space telescope 
( Appourc haux et al. 2008; Michel et al. 2008). The rich oscilla- 
tion spectrum of HD 49933 has recently been subject to some 
discu ssion, in particular in context with its funda mental param- 
eters (iBruntt etal .1120081; iRvabchikova et al.1 120091; Bruntt et al. 
2009, ac cepted) and the ability to model the oscillation fre- 
quencies (Kallinger et al. 2008a), underlining the importance of 
studying stellar oscillations in metal-poor objects with well con- 
strained parameters. Further recent detections of solar-like os- 
cillations in metal-poor stars i nclude the CoRoT timeseries of 
HD 181906 ([Fe/H] = -0.1) dMichel et al.ll2008l) and Doppler 
measurements of HD 203608 dMosser et all l2008l) ([Fe/H] = 
-0.65). 

The detection of p-modes and the measurement of the large 
and small frequency spacings in the high-overtone eigenspec- 



trum of 85 Peg A would put important constraints on its mass 
and age. If the star were as old as 12 Gyr or more, this would 
also place a meaningful lower limit on th e age of the Galaxy. For 
the se reasons, even before the work of Fernan des et al.l ([2002) 
and lD' Antona et al.l d2005). 85 Peg was selected as a primary sci- 
ence target for the MOST space mission, to search for p-modes 
in the metal-poor sub-dwarf component of the binary. The con- 
tinuous time coverage, high sampling rate and high photometric 
precision of the MOST photometry also make the data suited to 
search for other variability in the binary system. 

2. MOST observations & data reduction 

The MOST (Microvariability & Oscillations of STars) space 
mission is the first satellite desig ned to specifically conduct as- 
terose ismic observations of stars (Wal ker et al.l l2003; Mat thews! 
120071) . Equipped with a 15 cm Rumak-Maksutov telescope feed- 
ing a CCD photometer through a custom broadband optical filter 
(350-700 nm), MOST can collect rapid high-precision photom- 
etry of bright stars with nearly uninterrupted coverage for up to 
2 months. 

Depending on the brightness of the target of interest, MOST 
uses two different observing modes. In the primary mode, the 
target starlight enters a field stop about 1 arcmin in diameter and 
illuminates an extended fixed image of the telescope pupil on the 
science CCD detector using a Fabry microlens. For fainter stars, 
direct photometry of stars is obtained corresponding to tradi- 
tional CCD photometry with a FWHM of about 2-3 pixels. With 
a pixel scale of ~ 3 arcsec, the binary nature of 85 Peg remains 
unresolved (the integrated signal, however, being dominated by 
85 Peg A which is about 16x brighter than its companion based 
on the observed visual magnitude difference). 

MOST observed 85 Peg in 2005 and 2007. In 2005, the vi- 
sual brightness limit for the primary observing mode was V ~ 6, 
making 85 Peg (V = 5.75) the faintest star ever observed us- 
ing Fabry imaging. In early 2006, the satellite tracking CCD 
stopped functioning, most likely due to a strong particle hit. 
Consequently, guide star observations had to be relocated to 
the science CCD, limiting the exposure times for the program 
stars and constraining the brightness limit for Fabry Imaging. 
Therefore, 85 Peg was observed using direct imaging in the 
follow-up run in 2007. 

For both observation modes, data reduction p ipelines have 
been developed and are full y described by Reege net al.l 
(2006) for Fabry imag ing and iRowe et al. I (I2006F as well as 
Huber & Reegen (2008) for direct imaging. The main issue in 
MOST data reduction is the removal of stray light influences 
from the data. In all programs, the concept of decorrelation is 
used which is based on the comparison of intensities originat- 
ing from the star with intensities of background readings. The 
85 Peg datasets were re duce d according to the routin es described 
bv lReegen etal] (120061) and lHuber & Reegenl d2008l) . 

85 Peg was observed by MOST in 2005 for 25.6 days with 
almost no interruptions. The exposure time was 50 seconds with 
a sampling interval of 55 seconds. After outlier corrections in 
the reduction pipeline the final light curve consisted of 34,323 
measurements in total. The reduced light curve shows very good 
quality with point-to-point scatter of 0.48 mmag, making it well 
suited to search for p-mode oscillations. In 2007, MOST re- 
turned to 85 Peg for another run spanning a total of 25.5 days, 
this time with observations being performed in direct imag- 
ing mode with considerably shorter exposure times and images 
mostly being obtained every fifth MOST orbit. Due to the more 
difficult thermal conditions caused by target switching, the data 
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Table 2. Observation log of the MOST 85 Peg runs in 2005 and 
2007. Values in brackets for the 2007 run correspond to the light 
curve with 4 datapoint bins (in order to make the scatter compa- 
rable to the 2005 run). 



observation run 


2005 


2007 


start date [HJD-2451545] 


2085.43 


2830.81 


end date [HJD-2451545] 


2111.01 


2842.50 


total obs. time [d] 


25.58 


11.69 


number of exposures 


34,323 


5 738 


duty cycle [%] 


85 


14 


exposure time [sec] 


50 


13.25 [53] 


sampling time [sec] 


55 


25 [100] 


point-to-point scatter [mmag] 


0.48 


0.83 [0.51] 



are influenced by strong instrumental trends, which in particular 
affected the beginning and end of the run. Therefore, only the 
central part of the light curve spanning over about 12 days was 
used for the analysis. The data also show a slightly higher scat- 
ter than the earlier run, with an increased point-to-point spread 
of the light curve binned to about the same exposure time as the 
2005 data. An observation log for both runs is given in Table [2] 

3. The search for p-modes 

3.1. Frequency analysis 

Due to the data quality difference, the search for p-modes was 
only performed in the 2005 data, using Fourier analysis and least 
squares fitting. The sig nificances of identified peaks were es- 
timated with SigSpec (lReegenll2007l) which, based on a false- 
alarm probability, incorporates frequency, amplitude and phase 
information into the calculation. Significant frequencies are 
identified through consecutive prewhitening and least squares 
fits. The analysis was performed up to the Nyquist frequency 
of the data set (v Nyq - 9 mHz 780d~ 1 ). 

Figure [1] shows the Fourier amplitude spectrum in the fre- 
quency region where p-modes have been predicted. The ampli- 
tudes have been normalized by N/2 where N is the number of 
datapoints in the light curve. The largest peaks in the amplitude 
spectrum are residuals of the background variation at the orbital 
frequency of the MOST satellite (164 //Hz = 14.19d _1 ) and its 
harmonics, marked in Figure Q] by vertically dotted lines. Note 
that even the largest of these has an amplitude of only about 
35 ppm (0.035 mmag), and above a frequency of 2.5 mHz, there 
is only one with an amplitude as high as 15 ppm. There is also 
some power in sidelobes of the harmonics of the satellite orbital 
frequency, spaced by 1 and 2d -1 , due to the daily modulation of 
scattered earthshine. 

Aside from these instrumental frequencies, there is no appar- 
ent power excess above the noise level visible in the spectrum. 
The SigSpec analysis shows that, above 0.57 mHz (50 d~'), there 
are only two peaks (not associated with MOST orbital harmon- 
ics or 1 d -1 sidelobes) with a significance greater than 4 (corre- 
sponding to an amplitude S/N of about 3.5). One of these peaks 
can be seen in Figure[T]near 2.55 mHz. 

3.2. Estimation of the detection limit 

The MOST observations constitute the first large-scale dedicated 
observing effort to detect p-mode oscillations in the 85 Peg sys- 
tem. Despite the null-result, the data should allow us to draw 
valuable conclusions for future efforts by setting firm limits on 
the possible pulsation properties of 85 Peg A. 




Frequency [mHz] 



Fig. 1. Fourier amplitude spectrum of the 2005 MOST light 
curve of 85 Peg in the frequency region where p-modes have 
been predicted. Vertical dotted lines correspond to harmonics 
of the orbital frequency of the MOST satellite. The horizon- 
tal dashed line shows the 3.5 cr detection limit (corresponding 
roughly to a spectral significance of 4) for the noise level in this 
frequency region. The insert displays the window function of the 
light curve. 



Considering the average noise level of 4 ppm in the ampli- 
tude spectra at a frequency range where the p-modes are pre- 
dicted by theory, we can set an upper limit of pulsation in the 
MOST observations in Fourier domain to about 14-18 ppm, as- 
suming a 3.5-4.5 cr detection limit. Simulations of 11 coherent 
sinusoidal variations inserted into the MOST photometry spaced 
by Av = 1 70/iHz (a value lying in the range of the models cal- 
culated bv lD'Antona et alj 12 005)) and centered around 3.4 mHz 
(using the relation of Kjeldsen & Bedding; ( 1 1995b for the es- 
timated mass, radius and effective temperature of 85 Peg A) 
showed that peaks with amplitudes of 15 ppm are securely recov- 
ered using the same approach as in the real data using SigSpec 
Note that this limit is only valid under the assumption that the 
mode lifetimes are longer than the observation timebase and 
therefore can be modelled as coherent sinusoidal variations. 

Going one step further, we have also included the nature of 
the pulsation in the simulations. It is widely accepted that p- 
mode oscillations in sun-like stars are driven by turbulent mo- 
tions in the convective layers in the stars, causing a stochastic 
excitation and subsequent damping of the pulsation modes rather 
than a coherent oscillation. It has been shown that the lifetimes 
of the oscillations for sun-like stars are very sh ort, i.e. on the 
order of a couple of days (e.g. for cCenA, see iBedding et al.l 
12004 . while both short as well as substantially longer lifetimes 
(i.e. on the order of we eks or months) have been detected in 
giant and subgiant sta rs dStello et al.ll2006t ICarrieret all 120071: 
iKallinger et al.ll2008bl) . While the situation of subdwarf stars re- 
mains undetermined, it seems appropriate from these results that 
with a dataset as long as the MOST observations the effects of 
stochastically driven pulsation should not be neglected. 

The main effect when analyzing such signals in the Fourier 
domain is that the oscillation will not be described purely by the 
window function of the dataset, but rather by a series of peaks 
scattered around the true frequency value. The shape of this dis- 
tribution is commonly described by a Lorentzian profile with the 
FWHM of the profile being correlated to the mode lifetime and 
the height, in combination with the width, being r elated to the 
total energy of the pulsation mode (Houdek 2006). Hence, the 
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resulting amplitude of a mode in a Fourier spectrum is not only 
simply dependent on the maximum displacement in the time se- 
ries, but also on the amount of dampin g and the resulting mean 
variation in the data set. As given by Bau din et al.l d2005l) and 
Barb an et al.l d2007l) . the mean mode amplitude A can be ex- 
pressed as 



with T denoting the effective length of the dataset, r is the 
mode lifetime, H the height of the Lorentzian profile in power 
(ppm 2 ) and *K the height of the Lorentzian profile in power den- 
sity (ppm 2 Hz _I fl The effective dataset length T is evaluated as 
the reciprocal value of the area under the spectral window (in 
power), and in the case of the MOST 85 Peg observations gives 
T = 1.8895 Msec = 2 1.869 d. 

To investigate limits on the mean mode amplitude for 
85 Peg A, we simulated solar-lik e oscillations implementing the 
concept by Chaplin et al.l d 1997b adapted to the sampling of the 
MOST 85 Peg 2005 time series. Using the same input frequen- 
cies as in the simulations containing coherent signal as described 
above, about 1000 datasets were generated for mode lifetimes 
ranging from 1-26 d in steps of 1 d and for mean amplitudes 
from 1-40 ppm in steps of 1 ppm. Note that in order to keep 
the computation time reasonable and to solely study the ef- 
fects of different amplitudes and mode lifetimes, each simulation 
was computed as a single realization (i.e. the stochastic excita- 
tion/damping pattern was the same for each simulated data set). 
The mean mode amplitude (which was set to the same value for 
all frequencies) was simulated by setting the time-domain stan- 
dard deviation of a single realization to a fixed value A rms . Note 
that this value is not the same as the maximum peak height in the 
resulting Fourier spectrum, which is a priori unknown. In the fol- 
lowing, the standard deviation of the variation in the time domain 
will be denoted RMS amplitude, while the resulting peak height 
in the Fourier spectrum will be referred to as Fourier amplitude. 
For each simulation, the mean S/N ratio of the six highest peaks 
was calculated (i.e., we assumed 6 out of 1 1 frequencies to be 
a successful detection of p-mode pulsation). The results of the 
simulations are shown in Figure|2] 

Not surprisingly, the simulations show that oscillations with 
larger amplitudes and longer mode lifetimes would be easier to 
detect. For example, an RMS amplitude of 16 ppm of the os- 
cillations would remain undetected if the mode lifetimes are 
very short, even though Fourier domain noise would in princi- 
ple allow a 4 <x detection. Only for lifetimes getting closer to the 
length of the dataset are the expected S/N ratios as known for 
a coherent signal fulfilled (e.g., a coherent sinusoidal variation 
with a RMS amplitude of 15 ppm produces a peak amplitude of 
about 20 ppm, corresponding to a S/N~5 in the 85 Peg data). 
The positive influence of longer lifetimes, however, decreases 
rapidly as the lifetime approaches the length of the dataset, 
reaching a detection level of about 10 ppm at t > 10 d. Note 
that in Equation [TJ H is only representative of the power excess 
height for large values of T/t, i.e. when the Lorentzian profile 
is well resolved. This explains why a signal formally described 
by, e.g., H = 18 ppm can remain undetected for long lifetimes 
(i.e. small T/t) in the MOST photometry simulations. With this 
in mind, the following conclusions on the oscillation properties 
of 85 Peg can be drawn: 1) Independent of the amplitude or the 

1 Note that this relation is only valid when the amplitude spectrum 
is normalized by N/2, where N is the number of datapoints in the time 
series. 




5 10 15 20 25 

Mode lifetime [d] 



Fig. 2. Signal-to-Noise ratio of the six highest peaks of simulated 
solar-like pulsation using the MOST 85 Peg data as a function of 
mode lifetime and RMS amplitude in the time domain. Points of 
constant Lorentzian profile height H are shown as dashed green 
lines (from the top right to the bottom left: 50 ppm, 30 ppm, 
15 ppm and 5 ppm). The thick solid red line shows a power law 
fit to bins of data points with S/N=3.5-4.5. 



mode lifetime of the oscillations, the height of the Lorentzian 
profiles must be lower than about 20 ppm to remain undetected 
in the MOST data; 2) If the mode lifetimes in 85 Peg are very 
short (~ days), the RMS amplitude of the signal must lie below 
a level of about 12-18 ppm. If the lifetimes are long (~ weeks), 
this limit lies at approximately 10 ppm. 

In order to predict upper limits for the height of the 
Lorentzian profiles for 85 Peg A which are independent of the 
observation timebase, we converted the values derived from the 
simulations into power density. The top panel of Figure|4]shows, 
for each mode lifetime bin of the simulations, data points with 
S/N=3.5^4.5 and the corresponding power law fit in power den- 
sity. Combined with the Fourier amplitude limit of ~ 16 ppm set 
by MOST, this provides firm limits on the oscillation properties 
of 85 Peg A and could be used as a prior information to plan and 
analyze further observations with the goal of detecting solar-like 
oscillations in this star. 

3.3. Comparison with theoretical values 

In order to compare our upper limits on the oscil lation ampli- 
tudes with predictions, we follow the method by Miche fet al.l 
d2009h to calculate bolometric amplitudes, which are indepen- 
dent of the instrument being used. Using the transmission curve 
of the MOST filter and an effective temperature of 5600 K for 
85 Peg A, we calculate a response function for radial modes o f 
Ro = 4.787. Following the convention by lMichel et all (|2009), 
we then convert our upper RMS amplitude limits as derived in 
the previous section into bolometric amplitudes per radial mode 
as 

4 (6I\ 

Aboi,i=o = ^-lyJW- (2) 

where the last term denotes the observed RMS intensity fluctua- 
tions, which we identify as our definition of the RMS amplitude. 
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Fig. 3. Upper limits on oscillation properties of 85 Peg A as a 
function of mode lifetime, based on simulations using MOST 
photometry. Top panel: Upper limit on the height of the 
Lorentzian profile in power density. Bottom panel: Upper limit 
on the bolometric amplitude per radial mode. The open square 
shows the expected value based on theoretical scaling relations 
with 3 <x error bars. 



Theoretical values calibrated to observations indicate 
that bolometric photometri c amplitudes scale as (L/M) a T~^ 
dKieldserT & Beddingl 11995), with numerical simulations sug- 
gesting a value of a ~ 0.7 ( Samadi et al.ll2Q07h . U sing the solar 
referen ce value of 2.53 ±0.11 ppm for Aboi,i=o by iMichel et al.l 
(2009), the mass, luminosity and effective temperature (with a 
more conservative error of 100 K) of 85 Peg A yield an expected 
theoretical amplitude of 2.2+0.2 ppm. As an estimate for the ex- 
pected mode lifetime we use the recently suggested scaling rela- 
tion bv lChaplin etafl (|2009) proposing t qc T i, and with a solar 
reference value of t = 3.2 + 0.2 d dChaplin etal.1120051) arrive 
at t — 3.6 + 0.3 d for the expected mode lifetime of 85 Peg A. 

The bottom panel of Figure [3] shows the upper limits for 
85 Peg A as derived from the simulations converted to bolo- 
metric amplitudes per radial mode, compared to the theoreti- 
cal values. Clearly, the low value for the theoretical amplitude 
would in any case remain undetected in the MOST 85 Peg data. 
This is also true for the mean velocity amplitu de of 9 cms -1 
estimated by Houdek (cf. iD'AntonaetalJ^OOSl) for a model 
of 85 P eg A, which, using the relation by Kiel dsen & Bed ding 
( 1995), would translate in a bolometric luminosity amplitude of 
~ 1 ppm. Considering these estimated values, we are not able to 
draw any conclusions about the mode lifetimes for this star and 
can conclude that our results as upper limits are in agreement 
with predictions. 



3.4. Solar-type oscillations or wishful thinking? 

In a further investigation of the 2005 data, a relaxed spectral sig- 
nificance limit of 3 (corresponding to a S/N of about 3) was used 
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Fig. 4. Model fitting results as a function of mass using frequen- 
cies extracted from randomly generated datasets with a limit of 
S/N > 3. Line labels indicate the number of frequencies out of 
the total number of detected frequencies in a given dataset used 
for the fit. The fitting result using the real MOST dataset is also 
shown. 



to search for p-modes in the 85 Peg data. In this case, SigSpec 
identified several peaks which were mostly located in two dis- 
tinct frequency regions, one reaching from 1-3 mHz, the other 
from about 5-7 mHz. Eliminating all signals which are due to 
known instrumental artifacts and considering only values in the 
frequency range where p-modes are predicted, these frequency 
candi dates were fitted to a dense grid of 300000 stellar mod- 
els dGuenfher et al.|[T992T) . The quality of the fit was evaluated 
by a x 1 test, with a value of x 2 ^ 1 meaning that, on average, 
the model frequencies that are used in the fit are within the un- 
certainties of the observed frequencies (see Guenther & Brown 
2004 for details). This attempt resulted in a model fitting to about 
65% of the observed frequencies with ax 2 < 1. This model is in 
agreement with the previously estimated position of 85 Peg A in 
the HR-Diagram, with a large frequency separation of 172yuHz. 

To test the credibility of this model fit, we generated 8 simu- 
lated datasets, 4 using gaussian distributed random numbers, the 
other half by randomly shuffling the 2005 light curve. In both 
cases the same sampling, frequency analysis and model fitting 
procedure as for the real light curve was applied. The results of 
this test are shown in Figure [4] with the fit quality^ 2 given as a 
function of mass. As the Figure shows, there is indeed one ran- 
dom dataset where frequencies could be extracted that also fit 
with ax 2 < 1, indicating a similar fit quality as the best model 
derived for the real data. 

The result of this test underlines that when searching for 
low-amplitude p-mode oscillations, it is extremely important to 
discern real evidence from wishf ul thinking. In a recent anal- 
ysis of MOST data on Procyon, iGuenther et all d2008l) have 
demonstrated the risks when searching for comb-like structures 
in power spectra of sun-like stars using a priori assumptions. The 
test presented here is yet another example for this, showing that 
a S/N limit of 3 in amplitude (hence S/N > 9 in power, which is 
more widely used) is not sufficient to distinguish real signal from 
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peaks generated from random noise, and can lead to misleading 
results. 



4. Long periodic variability 

Although designed as a non-differential photometer, the MOST 
space telescope has proven to be remarkably stable over long 
timescales, making it also suited to search for long-periodic vari- 
ability. Instrumental trends, if present, can in many cases be 
identified, reconstructed and corrected by analyzing the satel- 
lites telemetry data such as board or preamplifier temperature. 

The 2005 light curve of 85 Peg shows a dominant long-term 
instrumental trend which is confirmed by inspection of both the 
board and preamplifier temperature data. All three datasets are 
displayed in panels (a), (b) and (c) of Figure[5] In a first attempt 
to correct the trend seen in panel (a), we fitted a polynomial func- 
tion to the data (dashed line). The residuals of this fit, low-pass 
filtered to remove residual variability with periods of one day 
and averaged into bins of length 101 min (about one MOST or- 
bital period), show two cycles of a variability with a period of 
about lid and an amplitude of 0.3 mmag (panel (d) of Figure|5]l. 
The 2007 data, which showed a less complex instrumental trend 
corrected by a simple linear fit, seems to confirm this variabil- 
ity detected in 2005, showing one full cycle of a light variation 
with a comparable period and amplitude (Figure [5] panel (f)). 
Combining the 2005 and 2007 data, a Fourier analysis yields a 
best fit period of 1 1 .6 d for the variation. 

In a later stage of the analysis, an attempt was made to cor- 
rect the 2005 data trend by directly using and fitting the vari- 
ations seen in the satellite telemetry dat£0. For this purpose, 
the telemetry data were binned to two orbital period means and 
heavily smoothed, so only long-periodic features remain present. 
These smoothed versions (shown as solid lines in panels (b) and 
(c) of Figure[5]l were then correlated to the light curve in two seg- 
ments of the run, and using the comparison set showing the high- 
est correlation coefficient (in both cases the preamplifier temper- 
ature, see panel (b) of Figure |5]l all correlations were removed 
using linear regressions. The resulting fit which was subtracted 
is shown as a solid line in panel (a) of Figure [5] The resulting 
light curve (panel (e) in Figure [5]l shows a much more complex 
behaviour, with the first cycle of the lid variability almost en- 
tirely removed. 

The considerable discrepancy between the results of the two 
background removal methods is obviously alarming. However, 
as the variability of the first reduction method was confirmed in 
the 2007 dataset, we feel that we have observed the ~ lid pe- 
riodicity in both epochs of data. We stress that without follow- 
up observations to confirm the observed variability we cannot 
rule out the possibility that the correlation between the 2005 and 
2007 epochs of data is purely coincidental, and that the variabil- 
ity changed considerably between 2005 and 2007 as illustrated 
with our second reduction method. It must be noted that the us- 
age of the decorrelation technique in this case assumes a lin- 
ear relationship of temperature and intensity variability, which 
might be an oversimplification of the problem. Indeed recent 
results on the low frequency variability of the MOST data on 
Procyon has shown evidence that telemetry data decorrelation 
seems to overcorrect existing intrinsic variability (Huber et al., 
in preparation). Thus we feel that the most likely possibility is 
that we have observed a common ~ lid variability in our two 



2 Note that the brightest star simultaneously observed in the 2005 
run is ~ 3 mag fainter and hence no comparison stars can be used for 
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Fig. 5. Panel (a): Original 2005 85 Peg light curve (every tenth 
data point is shown) together with different fits used to correct 
the instrumental trend (see text). Panel (b): Satellite preampli- 
fier temperature curve (every tenth data point is shown). The 
solid line shows a smoothed version. Panel (c): Same as panel 
(b) but for satellite board temperature. Panel (d): Binned and 
low-pass filtered 2005 85 Peg light curve derived by subtracting 
a polynomial fit (dashed line in panel (a)). Panel (e): Same as 
in panel (d), but derived by decorrelating variations found in the 
satellite telemetry data (solid line in Figure panel (a)). Panel (f): 
Corrected 2007 85 Peg light curve. Note that the y-axis scale for 
panels (d), (e) and (f) is the same, but differs from the scale in 
panel (a). 



epochs of data. In the following sections we present attempts 
to interpret the variability assuming that the 2005 light curve in 
panel (d) of Figure 5 is indeed intrinsic. 
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The variability could originate in either of the two stars in the 
binary, since the MOST photometry does not resolve the individ- 
ual components. If it is due to star A, then the amplitude of the 
variability is 0.3 mmag. If it is due to star B, then the amplitude 
we observe has been diluted significantly by the light of star A. 
The intrinsic amplitude would be approximately 5 mmag, based 
on the visual magnitude difference of 85 Peg A and 85 Peg B. 

We have considered two possible explanations of the period- 
icity: (1) light effects due to the possible binarity of component 
B, and (2) rotational modulation of the light of either star A or B 
due to surface spots. 

4.1. Binarity of 85 Peg B 

85 Peg B has long been suspected to be a multiple system due 
to the discrepancy between the dynamically determined masses 
and the magnitude difference between stars A and B. The shape 
of the light curves in panels (d) and (f) of Figure[5]is not consis- 
tent with a simple detached eclipsing binary, but it could be due 
to tidal distortions or backheating in a slightly eccen tric binary 
system with a period of 1 1 .6 d. Indeed. iGriffinl (l2004i) points out 
that the period of a putative 85 Peg Bb of less than a few weeks 
would be possible. 

Based on the location of 85 Peg B i n the HR diagram and 
by applying a mass-luminosity relation, ID'Antona et al.l d2005l) 
suggested the luminosities of a possible secondary system to be 
Lbo ~ 0.064 Lq and Lgb ~ 0.0085 Lq (with the luminosity of the 
primary La ~ 0.617 Lq). The masses of the components of the 
possible 85 Peg Bab binary are estimated to be M Ba ~ 0.5 M 
and M Bb ~ 0.2 M G dBach et alj 120071: ID' Antona et alj 120051) . 
Therefore, an orbital period P = 1 1 .6 d means a semi-major axis 
a ~ 0.09 AU (about 20 R e ). If the stars in 85 Peg B are on the 
main sequence, then this separation would be about 20x the ex- 
pected diameter of star Ba. 

The amplitude of the light variation is about 5 mmag (allow- 
ing for dilution by star A), which means the 85 Peg Bab binary 
must produce a peak-to-peak variation of about 1%. If the light 
variation is primarily a geometrical effect of tidal distortion, the 
light curve would be a double wave (with two maxima and two 
minima) during each binary orbit. The latter would mean the bi- 
nary period would be 2 x 11.6 = 23.2 d. However, for such a 
period, the semi-major axis of the orbit would be a ~ 0.14 AU, 
or about 30x the diameter of star Ba. At this separation, the geo- 
metrical distortion of star Ba by tides would be far below 1 % and 
could not account for the full amplitude of the observed varia- 
tion. 

The bolometric flux on the surface of star Ba due to star 
Bb (L Bb ~ 0.0085 Lq at a distance of about 0.1 AU) would be 
~ 1 .5 x 10 4 W/m 2 . The surface flux of star Ba (L Ba ~ 0.064 Lq; 
R ~ 0.6 R Q ) would be roughly 1 .4x 10 8 W/m 2 . Therefore the flux 
ratio due to backheating by a putative star Bb is only about 10" 4 , 
and the expected peak-to-peak amplitude would be in the order 
of 120ppm. This is two orders of magnitude smaller than the ob- 
served peak variation in brightness, if that modulation originates 
in component B of the 85 Peg system. 

Is it possible that the brightness ratio of stars A and B is not 
~16 due to the custom design of the MOST bandpass, and that 
variability due to star B would not be diluted as much as we have 
assumed? We tested this possibility by generating stellar models 
across a range of masses, for ages of 8 Gyr, and determining their 
flux ratios measured through the MOST filter (Table[3]). 

To change the flux ratio measured by MOST significantly, 
the temperature of star A would have to be considerably cooler 
than measured for this relatively bright star. Even if the the ratio 



Table 3. Model parameters for several possible configurations 
of the 85 Peg system and flux ratios as measured through the 
MOST custom filter. The deviations i n tem perature to the ob- 
served value given in ID'Antona et all d2005l) are given in terms 
of the standard deviation <x. A Q b s indicates the observed ampli- 
tude scaled with the corresponding flux ratio. The expected am- 
plitude due to backheating of star Bb on star Ba is about 0.1- 
0.2 mmag. 
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is lowered by a factor of 3 (last row of Table [3), the observed 
amplitude would still be (in the most optimistic case) by about a 
factor of 9 too low. Therefore, it is unlikely that illumination of 
a component Ba by a putative star Bb can easily account for the 
amplitude of the observed variation. 

4.2. Rotational modulation 

4.2.1. Determination of rotational velocity 

In order to interpret the variation in terms of rotational mod- 
ulation, a good knowledge of the projected rotational veloc- 
ity v sin i is essential. Search ing the literature we found mea- 
surements in the c a talogu es of Bernacc a & Perinottol dl970h and 
lUesugi & F ukuda (1970) giving vsinz values of 3 + 3 and < 
6kms ~ 1 , respectively, for 85 Peg A. M ore recentl y, [yVblff et al. 
(1985) measured v sin; 4km s _1 and lHald(ll994l) gives a value 
of 1 .8 + 0.6 km s _1 . The discrepancy between all these measure- 
ments is evident, and hence v sin i of 85 Peg A remains a rather 
poorly determined parameter. Note that due to the close separa- 
tion, no measurements for 85 PegB are available. 

With the goal to arrive at an independent vsirn measure- 
ment including statistically meaningful uncertainty estimates, 
we extracted high-resolution (R ~ 60000) spectra taken with 
the UVES spectrograph from the ESO archive which were pre- 
viously used to study star formation in the galaxy through abun- 
dance ratios (project number 076.B-0133). 

After the pipeline reduction and continuum normalization, 
synthetic spectra were computed using the p rogram synth3 
(Kochukhov 2007) in conjunction with VALD (|Piskunov et al.l 
Il995fc lKupkaetal.lll999HRvabchikovaetal.lll999l) . The fitting 
of the syntheses to the observed Ha line, which can be seen in 
Figure [6] resulted in logg = 4.5 and an effective temperature 
r eff = 5500K. 

Base d on this effective temp erature, MARCS model atmo- 
spheres dGustafsson et al.l 12008) were chosen for the subse- 
quent analysis. After the equivalent wid ths were determined us- 
ing a modified version (lTsvmballll996l) of the WIDTH9 code 
(lKurucz| [T993), a microturbulence velocity (Vmi C ) of 1.3kms _1 
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Fig. 6. Observed UVES spectrum (thick black line) and fitted 
synthetic spectrum with \ogg = 4.5 and T e ff = 5500 K (thin red 
line). 
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Fig. 7. The goodness-of-fit of observed to synthetic spectrum of 
one Fe line as a function of v sin i and v mac . Bright areas indicate 
a good fit. 



was determined by eliminating any abundance-equivalent width 
correlation. Specific Fe lines with highly accurate loggf (T. 
Ryabchikova, private communication) were used for the v sin i 
determination. An Fe abundance of — 0.35 ± 0.05 dex relative 
to the solar value (lAs plund et al.l l2005) and a radial velocity of 
34.6 km s" 1 were adopted to enable a good comparison of obser- 
vation and synthesis. Since the effects of macroturbulence veloc- 
ity (v mac ) and vsini on the goodness-of-fit in the spectrum are 
indistinguishable to the naked eye (see Figure [7]), we pursued to 
approximate the marginal distribution of v sin i. 

For this purpose, we calculated a 2D-grid of synthetic spec- 
tra, using the parameter values mentioned above, but with v sin i 
ranging from to 4kms~ I and v mac ranging from 1 to 4kms~ 1 
in steps of 0.05. Subsequently, we compared this grid of models 
to the four most clearly defined iron lines in the observed spec- 
trum. However, each of these lines was treated as an independent 
observation to overall create four probability distributions. The 
error bars of each observed data point, as delivered by the reduc- 
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Fig. 8. The marginal distribution of the v sin i parameter derived 
as described in 94.2.11 There is no evidence for a v sin i above 
2.5 km s _I according to the UVES data we used for our analysis. 



tion pipeline, were interpreted as lcr-confidence limits and the 
probability for each model value was calculated according to 



p(doj) = exp - 



(<io,i ~~ d\ 



M.i 



2ov 



(3) 



were c/o,i is the i-th observed data point, a?M,i is the correspond- 
ing model value, and <j\ the aforementioned lcr-uncertainty. The 
probability pm for each model then is the product of the proba- 
bilities for all data points. 

We calculated the marginal distributions of v sin i for each 
line by integrating out the v mac parameter. Since we used a dis- 
crete parameter space in our 2D-grid, we summed up the model 
probabilities over the whole parameter range of v mac for each 
specific value of v sin i, and normalized the distribution so that 
the sum of all resulting probabilities is unity. Finally, we com- 
bined the marginal distributions of all Fe lines by calculating the 
joint probabilities for all models and, again, normalized the re- 
sulting distribution. The final distribution is shown in Figure[8] It 
was subsequently used as a prioi0 for the v sin i parameter in the 
spot model analysis described in the next section. Judging from 
the distribution alone, we obtain v sin i = 1.40*J"|j! (with the un- 
certainties roughly corresponding to a 1 <x confidence limit). 



4.2.2. Spot model fitting 

To test the possibility of rotational modulation, we fit the 2005 
and 2007 light curves (panels (d) and (f) of Figure|5]l with simple 
sp ot models. We emp loy t he StarSpotz functionality discussed 
in ICroll et al.l (l2006al) and ICrollI d2006l) . whic h adopts circula r 
spots and applies analytic models described by iBuddingl i 1 9771) . 
StarSpotz has been used to measure the differential rotation pro- 
file of another MOST ta rget, the young active solar-type star 
k 1 Ceti dWalker et al.l 120071) and it was able to account for the 
photometric variations detected by MOST and constrain the dif- 
ferential rotation coefficient of the exoplanet-hosting star eEri 



3 In Bayesian probability analysis, the prior is a probability distribu- 
tion conveying already existing knowledge about a quantity. The infor- 
mation contained in the prior is used to constrain and assess the final 
probability of the outcome (see Bayesian inference textbooks for de- 
tails). 
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dCroll et al . 2006b). The reader is reminded at this point that in 
the following modeling efforts we explicitly assume that the vari- 
ation seen in panel (d) of Figure[5]is intrinsic, and disregard evi- 
dence for more complex variability as seen in panel (e) of Figure 

s 

A separate spot model was generated for components A and 
B of the 85 Peg system. Linear limb darkening coefficients for 
both stars were estimate d using the closest grid points to the 
p hysical param eters from D'Antona et al. (2005) in the catalog 
of lClaretl d2000t) ; these were u ~ 0.86 for 85 Peg A and u ~ 0.80 
for 85 Peg B. Note that considering the uncertainty of the ob- 
served variation and the fact that we are only intending to pro- 
vide a first approximation to model the observed variation, we 
refrain from using a more complex limb-darkening model in this 
application. We used a flux contrast between the spotted and 
unspotted photosphere of k w = 0.30 for 85 Peg A and k w = 0.35 
for 85 PegB. We are unable to independently determine k w , and 
small variations in k w will only serve to decrease or increase the 
size of the modeled spots. The value for 85 Peg A is appropriate 
for a G5IV star as it is c omparable to that o f the Sun and for the 
value chosen for a: 1 Ceti (Wal ker et al.l2007l) . We assumed a sim- 
ilar value for the K6-8V star 85 PegB. For the 85 Peg B model, 
the light curve was scaled to an amplitude of 5 mmag to allow for 
the dilution by star A. For our 85 Peg A spot-model fit we place 
a priori constraint on the resulting v sin i value (thus a joint prior 
information on the rotation period and the rotational inclination 
angle of the star) based on the values calculated in 94.2.11 For 
85 Peg B we do not use such a prior constraint since the v sin i 
value of this target is unknown. 

If the modulation is due to a single spot rotating on star A, the 
near-sinusoidal shape and low amplitude (0.3 mmag) of the light 
curves require this spot to be small (for the adopted flux contrast) 
and for it to be visible for the duration of the light curve. Our 
determination of v sin i rules out edge-on and near edge-on in- 
clinations for 85 Peg A; however, the uncertainty in v sin i is still 
significant and thus our spot models of this star could vary from 
near pole-on inclination angles to much more moderate values 
(i ~ 30°). Also, the observed modulation could be due to many 
spots or more complicated spot geometries. The spot model that 
follows should thus be in no way considered unique, but sim- 
ply an example of a possible one-spot spot model that could be 
causing the observed modulation. 

The best model for star A using the 2005 data features a spot 
of diameter 2° at a latitude of about 50°, while the 2007 data 
features a spot of diameter 2° at a latitude of about 46° relative 
to the equator. Our spot model has a rotational inclination an- 
gle of / ~ 21°. Note that since changing the inclination angle of 
85 Peg A does not significantly affect the x 2 of the spot model, 
our best fit model has an inclination angle that results in a v sin i 
similar to the mean of the distribution determined by our analy- 
sis in 94~2TT1 For the ra dius ofR A = 0.82 + 0.02 R Q - modeled by 
IrrAntona et al] (120051) for the range of possible a ges of 85 Peg A 
and cl ose to the value of 0.846 R Q estimated bv lFernandes et all 
(2002) - and a rotation period P — 1 1 .6 d, the equatorial rotation 
velocity of the star would be 3.6kms~'. The spot model inclina- 
tion results in a value of vsin/ of about 1.3 km/s (compared to 
our v sin i ~ 1 .40^ value determined in 94.2. 11 1. 

For star B, the intrinsic amplitude of the light variation 
should be closer to 5 mmag, requiring a larger spot if the contrast 
ratio between the spotted and unspotted photosphere is similar 
to that of other K6-8V stars. The best-fit model is presented in 
Figure [9] (note that the best-fit model for star A is very similar, 
except with a larger inclination angle and smaller spot diameter). 
In this case, the inclination of the star for the best fitted model 



is i si 11°. In 2005 the spot in the best-fit model has a diame- 
ter of 8° at a latitude of 27°. In 2007 the spot has a diameter of 
9° at a latitude of 20°. The equatorial rotation velocity of star B 
would be about 4.7 km/s and the resulting v sin i of the model is 
about 0.9 km s -1 . We could find no measurements of vsin/ for 
85 Peg B in the literature. 

5. Conclusions 

Nearly continuous photometry of the binary system 85 Peg ob- 
tained with the MOST space telescope in 2005 yields a null 
result for the detection of p-modes in the metal-poor subdwarf 
85 Peg A. Spanning over 25 days and reaching down to a noise 
level of about 4ppm in the frequency region of interest, the high 
precision data present the most complete photometry of the sys- 
tem collected so far and enable us to set an upper limit for the 
detection of such oscillations. Based on simulations of stochas- 
tically excited and damped p-mode signal, we conclude that the 
RMS amplitude of oscillations in 85 Peg A must lie below ~ 12- 
18 ppm if the mode lifetimes are short (i.e. on the order of days) 
and below ~ 10 ppm if the mode lifetimes are long (i.e. on the 
order of weeks). Based on these results, we provide upper limits 
for the Lorentzian profile height and the bolometric amplitude 
per radial mode as a function of mode lifetime and find these 
limits to be in agreement with theoretical scaling relations. Due 
to the fact that the expected theoretical amplitude is a factor 4— 
5 lower than the detection threshold, however, we are unable to 
comment further on the compa rison of theor e tical a nd observed 
amlitudes as recently done by iMichel et al.l d2008l) for CoRoT 
observations. 

We detect light variations in the 85 Peg system with a period 
of about lid and a peak amplitude of 0.3 mmag in our 2005 
epoch of data and confirm this variability in the 2007 epoch of 
data. However, as different methods to remove the instrumen- 
tal background in the 2005 data yield vastly different results, 
we cannot exclude the possibility that we do not know the ex- 
act shape and period of the 2005 variability. Nevertheless, as 
the 2005 and 2007 data yield similar results we believe it is 
likely that 85 Peg displayed consistent variability over these two 
epochs. 

The observed peak amplitude of the periodicity can in any 
case not be explained by effects of a possible undetected binary 
85 PegB. Ignoring correlations with telemetry data which exist 
but cannot be interpreted beyond doubt, we fitted simple models 
to the residual data under the assumption that the variation is 
caused by star spots modulated with the rotation period of one 
of the components. To facilitate the search for such models we 
analyzed high-resolution UVES spectra of 85 Peg found in the 
ESO archive and derived a new estimate of v sin i = 1.40^ [q for 
85 Peg A. While we are able to reproduce the light curves of both 
datasets well with single-spot models, it must be emphasized that 
the solutions presented here are by no means unique, and there 
is no evidence that rules out more complex models. 

Assuming that our simple models are realistic, this would 
indica te a rotation period of 1 1.6 d for either one of the compo- 
nents. |NQye£etal] ([1983) derived a rotation period for 85 Peg A 
of 23.8 days indirectly based on its chromospheric ac tivity index. 
The evolutionary models of iD'Antona et al.l (120051) allow ages 
from 8-14 Gyr, while those of iBach et al.l d2007l) result in an age 
of 8.4 + 0.5 Gyr for both stars in the 85 Peg system. Component 
B is expected to be coeval with star A and should share the 
same age. All of these age estimates suggest that the stars 
should have a slower rotation period than the observed variation 
shows. For example, the much younger active G5 dwarfs k 1 Ceti 



10 



D. Huber et al.: MOST photometry of 85 Pegasi 





Fig. 9. Spot models of the periodic light variation observed by MOST in the 85 Peg system if caused by 85 Peg B. Top: The best 
fitting one-spot configurations for 85 Peg B with the spot facing the viewer for 2005 (left) and 2007 (right) Bottom left: The top 
panel shows the observed (binned) 2005 light curve and the model fit for 85 Peg B, while the bottom panel displays the residuals to 
the fit. Bottom right: Same as the middle panel except for the 2007 light curve. 



(IWalker et al.ll2007h and HD 189733 dHenrv & Winnl2008h have 
rotation periods of 9 and 11.8 days, respectively - comparable 
to the 11.6 days rotation period suggested for 85 Peg despite the 
fact that the 85 Peg system is much older and should have a much 
longer rotation period. 

Th e spectroscopic binary solution summarized by Griffin 
(2004) yields an orbit inclination of | i \ ~ 50°. While we do not 
explicitly assume a spin-orbit coupling in the system to constrain 
the inclination in the spot models described above, it is worth 
noting that for 85 Peg A, a rotation period of twice the observed 
period (23.2 d) with such an inclination would result in a vsin/ 
of 1 .4 km s , which is also consistent with the values derived in 
94.2. H and might be a more realistic rotation period (suggesting 
that indeed a more complex spot model is required to explain 
the variation). Another, maybe more exotic explanation for the 
short period might be that the rotational modulation is produced 
in star B (the cooler and potentially more active component) and 
that its rotation is synchronized with the long suspected binary 
companion with an orbital period of 1 1 .6 d. 

Lastly, we strongly encourage follow-up observations to con- 
firm the 1 1 .6 d observed variability. 
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